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ABSTRACT 

We compare the metallicity distribution of Galactic Halo stars with 3D realizations of hierarchical 
galaxy formation. Outflows from dwarf galaxies enrich the intergalactic medium inhomogeneously, at a 
rate depending on the local galaxy density. Consequently, the first stars created in small early-forming 
galaxies are less metal-rich that the first stars formed in more massive galaxies which typically form 
later. As most halo stars are likely to originate in accreted dwarfs, while disk stars formed out of 
outflow-enriched gas, this scenario naturally generates a "metallicity floor" for old disk stars, which we 
find to be roughly coincident with the higher end of our predicted metallicity distribution of halo stars, 
in agreement with observations. The broad and centrally peaked distribution of halo star metallicities 
is well reproduced in our models, with a natural dispersion depending on the exact accretion history. 
Our modeling includes the important "baryonic stripping" effect of early outflows, which brush away the 
tenuously held gas in neighboring pre-virialized density perturbations. This stripping process does not 
significantly modify the predicted shape of the halo star metal distribution but inhibits star-formation 
and hence the number of accreted stars, helping to reconcile our model with the observed total Galactic 
halo luminosity and the lack of low- luminosity local dwarf galaxies relative to N-body predictions. 
Subject headings: Galaxy: halo - Galaxy: formation - intergalactic medium - galaxies: dwarf - 
cosmology: theory 



1. INTRODUCTION 

There is no completely successful model for the metal 
distribution of stars in our Galaxy. The simple 'one-zone' 
model for the disk, in which enrichment increases with 
time, is a poor match to the data, (van den Bergh 1962; 
Schmidt 1963) generating a broad and weakly peaked dis- 
tribution. Better agreement is found by extending the 
disk-formation time with a declining gas inflow rate (Lar- 
son 1974; Pagel 1989) but this cannot explain the sharp 
cut-off in the number of G-dwarfs with metallicities below 
[Fe/H] ~ — 1. In order to explain the full distribution, 
variable infall must be complemented by an ad hoc ini- 
tial production spike or "prompt initial enrichment" (e.g., 
Truran & Cameron 1971; Ostriker & Thuan 1975; Koppen 
& Arimoto 1990). 

In the Galactic Halo, a very different enrichment history 
is required. Here the metallicity distribution is broader 
with a peak shifted to a lower value of [Fe/H] ~ —1.5 
and is over-abundant in a elements compared to the disk, 
implicating predominantly supernova Type-II enrichment 
(e.g., Gilmore & Wyse 1998). No sharp cutoff is found but 
rather a long tail extending to low values of [Fe/H] = — 3 
and below (Ryan & Norris 1991). Unlike the disk, this 
metal distribution is generally understood as a sum of in- 
dividual dwarf galaxies accreted by the Galaxy. This idea 
was first suggested by Searle & Zinn (1978) and survives 
today in hierarchical models, such as the Cold Dark Mat- 
ter (CDM) model (e.g., Kauffman, White, & Guidcrdoni 
1993; Cole et al. 1994). Observationally, this scenario is 
strongly supported by tidally disturbed objects, in partic- 
ular the Sagittarius dwarf galaxy and its associated stream 
of Carbon stars (Ibata, Gilmore, & Irwin 1994; Ibata et 



al. 2000), by concentrations in the velocity and spatial dis- 
tribution of halo field stars (Preston, Beers, & Shectman, 
1994; Carney et al. 1996, Helmi & White 1998), and the 
large spread in the ages of both the halo field and halo 
cluster stars (Laird & Sneden 1996; Sarajedini, Chaboyer, 
& Demarque 1997). For a different point of view, see how- 
ever, Oey 2000. 

Despite the growing consensus that the Halo was assem- 
bled hierarchically, a realistic description of this process 
has not been developed. Argast et al. (2000) consider a 
3D model of interstellar medium evolution that roughly 
matches the observed properties of the Halo but is un- 
connected to the accretion of dwarf objects (see also Tsu- 
jimoto, Shigeyama, & Yoshii 1999). Cote et al. (2000) 
adopt an empirical approach, constructing the metal dis- 
tribution from the observed metallicity- luminosity relation 
of low-mass galaxies. This model, while roughly successful, 
is at odds with the observed paucity of local dwarf galax- 
ies (Klypin et al. 1999; Moore et al. 1999) and leaves open 
the issue of the origin of metals in the dwarf objects them- 
selves. Finally, Hernandez & Ferrara (2000) examined the 
lowest metallicity halo stars in terms of the mass- weighted 
average intergalactic medium (IGM) metallicity for each 
dwarf object formed, an approach which is most reliable 
for lower metallicity values. 

In this letter we offer a holistic 3D picture in which 
the physics behind the luminosity-metallicity relation, the 
cause of "prompt initial enrichment" of the Galactic disk, 
and the mechanism that prevents an excess of local group 
dwarfs all stem from the same process. Galaxy outflows 
enrich the IGM slowly and inhomogeneously over time, 
so that the first stars formed in small early galaxies are 
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more metal poor than the first stars formed in larger, later- 
forming objects. As prc-virialized overdensities that would 
otherwise form galaxies are especially susceptible to dis- 
ruption by winds from neighboring objects (Scannapieco, 
Ferrara, & Broadhurst 2000; Scannapieco & Broadhurst 
2001, hereafter SB; Scannapieco, Thacker & Davis 2000) 
galaxy outflows preferentially suppress the formation of 
late-forming dwarf galaxies, reducing the total number of 
local group dwarfs and Galactic Halo stars. 

The structure of this work is as follows. In §2 we de- 
scribe our semi-analytical simulations of galaxy formation 
and metal dispersal, and our criteria for selecting MW-like 
objects. In §3 we present the results of these simulations 
in relation to the Halo, with conclusions in §4. 

2. MODELS OF THE MILKY-WAY STELLAR DISTRIBUTION 
2.1. Simulation of Galaxy Formation 

Using the semi-analytical model described in detail in 
SB, we construct a (12 Mpc//is) 3 cubic comoving volume, 
where h is the Hubble constant in units of 100 km/s/Mpc, 
and track the collapse of dark mater halos, cooling of gas, 
and formation of galaxy outflows. A (512) 3 linear density 
field is constructed from a standard fit to the galaxy power 
spectrum and convolved with spherical "top-hat" filters 
corresponding to ten different mass-scales, spaced in equal 
logarithmic intervals from 1.1 x 10 8 M Q to 5.3 x 10 12 M Q 
and covering the interesting range from the most massive 
galaxies to objects that lie close to the lower limit set by 
photo- ionization and molecular cooling (SB). Using the 
elliptical model for collapse developed in Sheth, Mo, & 
Tormen (1999) we search through these fields, identifying 
the collapse redshifts and locations of dark-matter halos. 

Additional modeling is required to reproduce the forma- 
tion of galaxies. The two important issues in this case are 
the delay between halo collapse and the cooling of the viri- 
alizcd gas into stars, and the ejection of material into the 
IGM by the outflows generated in the initial star-formation 
episode. 

To account for the delay between collapse and star- 
formation, we use the standard inside-out collapse model, 
calibrated against numerical simulations (White & Frenk 
1991; Somerville 1997). The cooling gas is modeled as an 
initially isothermal sphere (p(r) cx r~ 2 ) at the virial tem- 
perature T and with a constant metallicity, Z calculated 
as described below. All the gas within a radius r coo i cools 
instantaneously according to a p 2 A(T, Z) radiation law, 
where A is the radiative cooling function, and the gas out- 
side this radius stays at the virial temperature of the halo, 
with r coo i moving outward with time. When the total mass 
contained within r coo i equals the object's baryonic mass, 
we consider a new galaxy to be formed. 

Once an object has collapsed and cooled into a galaxy, 
we assume that a fixed fraction of the stellar material e s f 
is converted into stars in an initial burst of star-formation 
and that an outflow develops powered by supernovae and 
stellar winds. We estimate that one massive star forms for 
every 100M Q of stars, with an associated energy output 
of 1 x 10 51 ergs in stellar winds during its lifetime and an 
additional 1 x 10 51 ergs when it explodes. A fraction e w i n( j 
of this energy is channeled into a galaxy outflow, ejecting 
gas into the IGM. 

The outflows are modeled as spherical shells that are 



driven by internal pressure and decelerated by gravita- 
tional breaking by the halo (see SB), both estimated in the 
thin shell approximation (Ostriker & McKee 1988). These 
shells expand into the Hubble flow, sweeping up the IGM 
and losing only a small fraction f m = 0.1 to the interior. 
Their evolution is then completely determined by the as- 
sumed luminosity history of the blast wave, and we follow 
Tegmark, Silk, & Evrard (1993) in considering only four 
contributions: energy injection from supernovae, heating 
from collisions between the shell and the IGM, cooling 
by Compton drag against the microwave background, and 
bremsstrahlung radiation and other two-body interactions. 

The metallicity of each collapsing object is volume aver- 
aged over the contributions of each of the expanding blast 
waves passing within the collapse radius r co ii, adopting 
an average yield from each supernovae of 2 M Q in metals 
(e.g., Nagataki & Sato 1997). Each object is assigned a 
mass in metals Mz, taken to be zero initially and mod- 
ified by each outflow passing within r co \\. For each such 
occurrence, Mz is updated to 
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where Vo VC riap is the volume of intersection between the 
outflow and the collapsing sphere. By dividing this mass 
by the total baryonic mass of the galaxy we can compute 
the initial metallicity of the object. 

Our outflow models incorporate the "baryonic strip- 
ping" mechanism identified in Scannapieco, Ferrara, & 
Broadhurst (2000), whereby outflows from early galaxies 
strip the gas out of nearby overdense regions that would 
have otherwise become relatively low-mass galaxies. Strip- 
ping occurs whenever a shock moves through the center of 
a pre-virialized halo with sufficient momentum to acceler- 
ate the gas to the escape velocity, v e , and thus we exclude 
all objects for which 



fM s v s > M c v e , 



(2) 



where / is the solid angle of the shell subtended by the col- 
lapsing halo, M s is the IGM mass swept up by the shock, 
v a the velocity of the shock, and M p is the baryonic mass 
of the pre-virialized halo. 

2.2. Identification of Milky- Way Type Galaxies and 
Construction of Stellar Distributions 

Having developed a model of galaxy formation, we now 
select MW-like galaxies from this simulation. Mass esti- 
mates of our Galaxy from line-of-sight satellite velocities 
place a lower limit of w 10 12 M Q (e.g., Kulessa & Lynden- 
Bell 1992; Binney & Merrifield 1998) within « 50kpc/h. 
For the elliptical collapse model employed in our calcu- 
lations the total mass is greater however, corresponding 
to the larger virial radius r 2 oo, where the average interior 
overdensity is 200. N-body comparisons between circu- 
lar velocity and total halo mass within r2oo suggest values 
~ 2.5 x 1O 12 M (e.g., Steinmetz & Navarro 1999). We 
adopt this approximately and select halos of 1.6 x 1O 12 M , 
the closest mass-scale in our simulations and a value simi- 
lar to 2.0 x 10 12 M Q used in Hernandez and Ferrara (2000). 

A second important property of the Milky Way is that 
it has not had a major merging event. The morphology- 
density relation (e.g., Abell 1958; Balland, Silk & Schaef- 
fer 1998), as well as numerical studies of galaxy formation 
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(e.g., Aguilar & White 1985; Hernquist 1993) suggest that 
the spherical components of galaxies form by major merg- 
ers. As the bulge of our Galaxy represents <J 10% of the 
total mass (e.g., Rohlfs & Kreitschmann 1988; Dehnen & 
Binney 1998) we therefore exclude all objects with more 
than one progenitor on the two mass-scales below that of 
the MW (4.8 x 10 n M Q and 1.4 x lO n Af ). 

The final constraint on the identification of MW candi- 
dates is the age of the assembled product, which should 
be relatively large to match the age of the MW disk. This 
is important as the IGM is formed over a protracted pe- 
riod, and thus a galaxy that formed much earlier or later 
is likely to have a different metallicity. Note that the time 
at which the collapsed gas within the 1.6 x 10 12 M Q object 
cools is associated with the formation of the Galactic disk 
only. A reliable measure of this age is given by the lower- 
limiting luminosity of white dwarfs in the disk, which cool 
with time (Schmidt 1959). Improvements in both obser- 
vations and our understanding of white dwarf structure 
have led to a consensus of approximately 10 Gyr for the 
age of the disk (e.g., Winget et al. 1987; Leggit, Ruiz, & 
Begeron 1998; Knox, Hawkins, & Hambly 1999). In our 
simulations we allow all objects that formed 10 ±3 Gyrs 
ago to be considered as MW candidates. 

Within each MW-like galaxy, we use the full merger his- 
tory to divide its stellar population into disk, bulge, and 
halo stars. The first generation of disk stars are tagged 
as those formed in the MW-like object itself, while the 
other two populations are formed within progenitor ob- 
jects. Whether progenitor stars are absorbed into the 
Galactic Halo or instead become part of the bulge depends 
on their initial trajectories. Here we apply the rough rule 
of thumb that the bulge was formed out of objects con- 
tained in the inner 10% of the Galactic dark-matter halo, 
while the stars in objects in the remainder of the collapsed 
region were accreted but not absorbed into the central 
bulge. As the orbits of these objects are much different 
than that of the main disk component, they are naturally 
associated with the Halo (e.g., Bekki & Chiba 2000). 

Note that the relatively early formation times of the 
MW candidates, corresponding to even earlier virialization 
times for their dark matter halos, mean that the majority 
of the the progenitor galaxies used to construct the halo 
distribution are formed at relatively high redshifts gtrsim 
2. These early times, coupled with a large minimum pro- 
genitor mass set by photo-ionization and molecular cool- 
ing concerns, insures that our semi-analytical treatment 
can be trusted to construct an accurate merging history 
for each candidate. 

Note also that our simulations only provide the metallic- 
ity of the first generation of stars to form in each object, 
and no subsequent star-formation is considered. Fortu- 
nately all accreted objects are sufficiently small (<; 5 x 
1O 1O M ) that second generation star-formation is quashed 
because much of the gas is expelled. This is consistent with 
observations, which suggest that only a small fraction of 
stars in the halo were enriched by Type IA supernovae 
(Gilmore and Wyse 1998). Hence the observed halo stellar 
metal distribution can be reliably compared to our simula- 
tions. Because our filtering quantizes the galaxy mass dis- 
tribution, we incorporate a small smoothing of the metal- 
licity of each galaxy, modeling it with a mean metallicity 



as calculated from Eq. (|]) and a dispersion of er([Fe/H]) 
= 0.3 independent the mass and mean metallicity of the 
dwarf (Cote et al. 2000), a procedure that introduces no 
significant bias and is more realistic than adopting a single 
metallicity for all the stars in a given object. 

3. RESULTS 

With these prescriptions we generate and identify model 
MW-like galaxies, which are massive, old, and have suf- 
fered at most a merger with less than 10% of their to- 
tal mass. We restrict our attention to the currently fa- 
vored ACDM cosmology and fix £7o = 0.35, f^A = 0.65, 
n b = 0.06, a s = 0.87, T = 0.18, and h = 0.65, where O , 
Q\, and Of, are the total matter, vacuum, and baryonic 
densities in units of the critical density, erg is the ampli- 
tude of mass fluctuations on the 8 h~ 1 Mpc scale, and T 
is the CDM "Shape Parameter." 

In our fiducial outflow model we adopt the same pa- 
rameters used in SB and Scannapieco, Thacker, & Davis 
(2000) and take e s f = 0.1 and e W i n d = 0.1. For comparison 
we also consider a "low-energy outflow" simulation with 
e s f = 0.1 and e w ; n d = 0.05, and a simulation in which we 
again set e s f = 0.1 and e w i n d = 0.1 but do not apply the 
baryonic stripping criteria, Eq. (||). In the fiducial and 
low-energy outflow models, roughly two thirds of the ob- 
jects in the 10 9 M Q to lO n M mass range are suppressed 
by stripping, helping to reconcile CDM predictions with 
the observed lack of visible local group dwarf galaxies. In 
all cases, the average initial metallicity of galaxies ranges 
from approximately [Fe/H] = -2.0 to [Fe/H] = -0.6, in- 
creasing roughly linearly with the log of the mass. Further 
details are presented in SB. 

The metallicity distribution of halo stars in the five MW- 
like galaxies identified are shown in Figure 1, along with 
the observations of MW Halo stars by Ryan and Norris 
(1991). As each simulation shares the same random seed, 
these model distributions are directly comparable across 
runs. In all cases, we see that the initial disk metallic- 
ity coincides with the high-metallicity end of the halo star 
metal distribution, mimicking the behavior of the MW. 
While the low-energy outflow and fiducial models produce 
similar results, reducing e w i n d = 0.05 slightly lowers Z val- 
ues, more in line with the observations but well within the 
uncertainties of our approach. Note that the metallici- 
ties in our models can be shifted linearly by altering the 
assumed stellar yields. 

Exclusion of baryonic stripping has a larger effect, pro- 
ducing a sharper peak and a factor of ~ 3-5 more halo 
stars. While the overall shape of the distribution is reason- 
able, this model suffers from the same excess of low-mass 
dwarf galaxies as the empirical model suggested in Cote et 
al. (2000). The models that include stripping via Eq. (2), 
on the other hand, can be more easily reconciled with the 
observed luminosity of the Halo: Ly w5x 10 7 Lq (Binney 
& Merrifield 1998) corresponding to a total stellar mass of 
~5x 10 8 Mq for a Scalo (1986) initial mass function and 
a ~ 10 10 year old stellar population. Note that as more 
winds are formed at later times, suppression is most severe 
in the later forming MW-like galaxies. 

Finally all distributions contain an excess of [Fe/H] 
~ —3 stars, corresponding to unpolluted (Pop III) galax- 
ies. This discrepancy most likely relates to our simple 
treatment of the star-formation within dwarfs, which fails 
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to account for what is undoubtedly a much altered pro- 
cess under primordial conditions (e.g. Haiman, Thoul, & 
Loeb 1996; Nakamura & Umemura 1998). Note that this 
discrepancy is qualitatively in agreement with Hernandez 
and Ferrara (2000), although, in our case the distribution 
of low-metallicity stars is somewhat complicated by our 
imposing at all redshifts the finite minimum mass scale 
set by photo-ionization. 

4. CONCLUSIONS 

We have considered the formation of the Galactic Halo 
in the context of hierarchical structure formation mod- 
els. While extensive theoretical and observational evidence 
suggests that the Halo was formed primarily from accreted 
dwarfs, a detailed understanding of the relationship be- 
tween this process and the metallicity distribution of the 
Halo has yet to be developed. 

Using a simple semi-analytic model, we have investi- 
gated the qualitative features that can be expected in a 
complete picture. Galaxy outflows enrich the IGM slowly 
over time, causing the stars in early-forming galaxies to be 
more metal poor than those in more massive later-forming 
objects. This results in a "metallicity floor" in the distri- 



bution of Galactic disk stars at a level roughly coincident 
with the high end of the distribution of halo stars, as ob- 
served. While previous models of hierarchical formation 
are at odds with the observed lack of local group dwarf 
galaxies, accounting for "baryonic stripping" allows us to 
construct realistic metallicity distributions while reducing 
the number of low-luminosity objects. 

Our models cannot claim to be complete, but rather 
point to the features that can be expected in a fuller hy- 
drodynamics! treatment. The sensitivity of our results to 
the strength of the outflow winds and the degree of bary- 
onic stripping suggests that a detailed understanding will 
require more careful numerical modeling of structure for- 
mation and outflows which must be consistent with the 
degree of halo velocity and metallicity substructure inher- 
ent to accretion. It is clear that only when we have un- 
derstood the enrichment history of the Universe and its 
relation to structure formation can we hope to develop a 
completely successful model for the metal distribution of 
stars in our Galaxy. 

We thank Andrea Ferrara and the referee, Xavier Her- 
nandez, for helpful comments and discussions. 
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Fig. 1. — Distribution of metallicities in the halos of MW-like objects. The solid lines are the results of the fiducial e w i n( j = 0.1 outflow 
model, the dotted lines are from the low-energy outflows simulation e w i n d = 0.05, and the dashed lines are the result of including halos which 
fail to pass the baryonic stripping criteria, Eq. (Eh, in the e w i n( j = 0.1 model. Each panel is labeled with the formation time, halo mass in 
units of IO^Mq, and initial metallicity of the disK component in the fiducial model, which is also labeled by the vertical lines. In the lower 
right panel we plot the metallicity distribution of the sample of MW Halo stars studied by Ryan and Norris (1991), arbitrarily normalized. 



